Universalities in vortex transport at the melting transition 
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We consider the jump in resistance at the mehing transition, which is experimentally observed to 
be constant, independent of magnetic field (vortex density). We present an explanation of this effect 
based on vortex cuttings, and universalities of the structure factor at the freezing transition (the 
Hansen- Ver let criterion) . 
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The subject of vortex dynamics in high Tc supercon- 
ductors is of topical interest ||l|,|[. In recent years, new 
aspects of the phase diagram of the vortex lattice have 
been elucidated through experiment. Perhaps one of 
the salient recent experimental findings is the first or- 
der melting transition of the vortex lattice, first observed 
in transport measurements O] and later confirmed in a 
number of equilibrium properties B. It is now well ac- 
cepted that in clean BiSrCuO and in untwined YBaCuO 
samples the melting transition of the vortex lattice for 
an applied magnetic field parallel to the c-axis is first or- 
der and that the superconducting coherence is simulta- 
neously lost in all directions at the melting temperature. 
A striking regularity is observed in the transport experi- 
ments: for different values of the applied magnetic field, 
which is proportional to the areal density of vortices, the 
in-plane resistivity Pab{T) decreases with decreasing tem- 
perature until it jumps (essentially) discontinuously to 
zero at the melting temperature Tm, with Pab{Tm) inde- 
pendent of field. Moreover, recent studies of the melting 
transition in YBaCuO show that the jump in the out-of- 
plane resistivity pc is also field independent. 

This universality in the jump in resistance at the melt- 
ing temperature constitutes, we argue, an evidence of 
general universalities in the structure factor S{q) of the 
liquid phase at the melting point as first pointed out by 
Hansen and Verlet for normal liquids M . 

Verlet's criterion for freezing is the counterpart of Lin- 
deman's criterion for melting. Verlet observed that the 
rescaled curves of S{q) are nearly identical for a variety 
of liquids that freeze to the same crystalline structure. 
The only scale is q™, the wave vector at the first peak in 
S{q), with S{qm) ^ 2.85. 

While Lindeman's criterion tells us a way of antici- 
pating melting by looking at mean square displacements 
of the atoms in the solid phase Q, Verlet's criterion 
provides a way of predicting freezing by looking at the 
structure factor in the liquid phase. Ramakrishnan and 
Yussouf iQ have offered a theoretical explanation of the 
criterion based on a phenomenological order parameter 
theory. 



More recently Lowen et al. S made the empirical ob- 
servation of a dynamical criterion for freezing of colloidal 
liquids. For these systems the ratio of the long-time and 
short-time diffusion coefficients of the colloidal particles 
Dl/Ds has a universal value close to 0.1 at the freezing 
temperature of a variety of systems. It has been shown 
that, for the colloidal liquids, universalities in S{q) at 
freezing imply universalities in the long time diffusion 
constant ||,|l^. 

In the present paper, we will assume the validity of 
Verlet's criterion for vortices. However, we wish to em- 
phasize that there are still important differences in the 
the diffusion coefficients of vortices and colloidal parti- 
cles which stem from the different topological properties 
of these systems. Vortices are extended objects that can 
entangle. In the entangled phase diffusion implies cutting 
of vortex lines, an effect which not present in suspended 
particles. For clean systems we show that vortex cutting, 
together with Verlet's criterion, implies a universality in 
the jump of both pab and pc- We illustrate this by solv- 
ing a simple model for vortex cutting, and show that the 
resistivity at the melting temperature is given by the nor- 
mal resistivity and the structure factor only, and thus the 
observed universality in the resistivity jump implies the 
validity of Verlet's criterion. 

We will first analyze the out of plane resistivity p^- 
The mechanisms of dissipation when the current is par- 
allel to the applied field originate in the entanglement of 
the field-induced vortices and the thermally generated 
vortex loops psl. In this entangled phase, the structure 
of vortex lines percolates through the sample in the di- 
rection perpendicular to the external field. The result is 
a "phantom mesh" of vortex lines that threads the sam- 
ple in all directions and where each vortex diffuses by 
cutting the others. The topology of the phantom mesh 
allows different segments of the field-induced vortices to 
diffuse independently. This implies that the vortex mo- 
tion is uncorrelated in the c-direction, as observed in the 
pseudo DC transformer experiments |14| . 

When a current flows along the c-axis, only those lines 
oriented in the direction perpendicular to the c-axis will 



contribute to the bulk resistivity. The resistivity can be 
written as 
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with np the density of percolating vortex lines, and Dp 
their diffusion coefficient. In the computation of Dp one 
needs to include the cutting of the percolating vortex 
lines with the field-induced vortices. 

We compute the diffusion by considering the random 
walk of the center of mass of a vortex line. We introduce a 
discretization of the random walk, with step length i, and 
we call Ti the time elapsed at the z-th step. Neglecting 
vortex cutting all time steps are equal, with Ti — tq, and 
the bare diffusion coefficient is given by Do — I'^/tq. This 
bare diffusion constant is related to the Bardeen-Stephen 
||12| , [ll| viscosity coefficient 770 = kpT/Do- 

Assume now that there is a fraction ric of steps, at 
which vortex cutting takes place, with a characteristic 
step time t^ ; the diffusion coefficient is then 



Dp = f/ [ricTy, + (1 - nc)To] 



The cutting frequency is 
1 1 



Tx To 



-Uy/keT 



(2) 



(3) 



with [/x an effective barrier for cutting iTS] . In the above 
equation we have assumed thermal activation for the vor- 
tex crossings. The Boltzman factor gives the relative 
probability of finding a vortex in a crossing configura- 
tion. 

li ksT < Ux , and assuming that Uc is of the order of 
the density of field-induced vortices n, the resistivity is 
given by 
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where ao 



~^/^ is the mean distance between field in- 
duced vortices along the c-direction and ^2 = Up is 
the mean distance between percolating vortices in the 
direction perpendicular to the field. In what follows, we 
argue that Ver let's criterion implies that pc is universal 
(in a weak sense) at the melting temperature. The ar- 
gument is two-fold: i-The universality of the structure 
factor S{q) indicates that the scaled correlation lengths 
are independent of density at Tm , and therefore the factor 
{o-o/^z) is density independent at freezing. Ji-The expo- 
nential factor in Eq. fl) represents an effective probabil- 
ity of cutting which will be also universal at melting due 
to Verlet' criterion. (See the solution of the simplified 
one-dimensional model below.) 

We conclude that the out-of-plane resistivity pc is in- 
dependent of the vortex density at the melting temper- 
ature Tm- This weak universality is a consequence of 



Verlet's criterion. The jump in the in-plane resistivity is 
also independent of the external field. The resistivity pat 
can be written as 
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with n = B/(f)Q the density of field induced vortices along 
the direction of the applied field B, D is the correspond- 
ing diffusion coefficient. The diffusion coefficient D can 
be calculated following the out-of-plane case, with an 
analogous results. We stress that the diffusion constant 
is reduced by a factor proportional to n due to vortex 
cutting, canceling the prefactor of n in the above equa- 
tion. 

The arguments presented above are valid for good 
quality single crystals for which vortex pinning is weak. 
Evidence from both experiments and numerical simula- 
tions have shown that intrinsic disorder destroys the first 
order transition. For example, upon irradiation, the mag- 
nitude of the jump decreases. In the cleanest YBCO 
samples, the jump in pab is about a fifth of the normal 
resistivity at the critical temperature. A smaller jump 
in the resistivity at the melting transition is an indica- 
tion of disorder. In our formulation, a large disorder can 
be incorporated in a modified particle-particle correla- 
tion function which will retain its universal properties as 
long as the transition remains first order. Consequently, 
to the extent that Verlet's criterion remains valid, the 
universality in the resistity jump is preserved. 

Now let us turn to a simplified one-dimensional model 
of interacting diffusing particles. The purpose of this dis- 
cussion is to justify using a probability of cutting which is 
universal at melting. We consider the following Fokker- 
Planck equation, describing diffusing interacting particles 
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with F, = -6u\f + Y.'^=^y'{xj-xo)\ + (1 - 
5i.o)V' {xi — XQ)^ V being an interaction potential between 
particle zero and the rest of the N particles. Note that 
there is no interaction between particle i and particle j 
for i ^ Q ^ j. Also, the diffusion constant is Dq for par- 
ticle 0, and Di -C -Do for z 7^ 0. We added a drag force / 
acting on particle zero, and computed the resulting mean 
drift velocity of particle zero {v) by solving the Focker- 
Planck equation imposing periodic boundary conditions 
||l6| on a length L. The result is {v) = {D/kBT)f, with 
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Note that, since particle zero diffuses faster tfian tfie 
other particles, we take the configuration of the rest of the 
particles as described by the equilibrium distribution in 
the absence of the force /. Also, if the potential is mono- 
tonically decreasing with distance with V{Q) » fcsT and 
the range of the interaction is much smaller than the 
inter-particle separation. 
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with A the range of the interaction. We can now iden- 
tify the integral in Eq. (|^) which involves the nega- 
tive exponential as a partition function Z , and defining 
Px = [exp-ViO)/kBT]/Z, we obtain 
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In this simplified model Px represent the probability of 
finding particle zero on top of particle i. This prob- 
ability can be extracted from the dimensionless pair 
correlation function g{r) or its Fourier transform S{q): 
P^=Jdiq/qm)[Siq/qra)-l]. 

For the vortex case the quantity nX of the above one- 
dimensional calculation is identified with n^^, ^ being 
the core radius. Finally, if in our equation for the re- 
sistivity pab = {(l)o/cYnD/kBT = (0o/c)^Px/(%?^) we 
replace the expression for the Bardeen-Stephen viscosity, 
770 = {h/ e)'^ {tt /2) / {pabN^.'^) with pabN the normal state 
resistivity , we obtain for the resistivity jump at Tm 
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where Px [S{q/ qm)] means that Px is a functional of 
S{q/qm) and is universal at melting. This equation is 
analogous to Eq.(^) and constitutes the main result of 
this work. 

The calculation of the probability Px of finding vor- 
tices in the cutting configuration requires the knowledge 
of the of the exact correlation functions g{r). An esti- 
mate can be extracted from the calculations of the "cage" 
model described in Ref . [Q . Within that simple an entan- 
glement length Z2 is defined as the distance along the field 
required for a vortex to diffuse a distance of the order of 
the vortex-vortex separation ao: Iz '^ g4>o/ikBTB), with 
g the tilt modulus, and Oq ^ 4>q/B. From the above 
discussion, we estimate 
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On the other hand, within the cage model, a melt- 
ing temperature can be extracted by applying Lin- 
demann's criteron on the solid phase [n|: /csTm = 
cj^£oirn±/mzy^^{(j}o/B)^^'^, with m± and rriz being re- 
spectively the in-plane effective mass and the out-of- 
plane effective mass, and eg = {rnz/m±)g a coupling con- 
stant giving the interaction energy per unit lenth for the 



vortex lines. Also, Cl is the Lindemann constant, which 
for vortices is of order 0.1-0.2 g]. Collecting this with 
Eq.(pll) we obtain 



Px - ci(mi/m,)i/2. 



(12) 



If we assume that the anisotropy on the effective mass 
is the same as that of the normal resistivities, we obtain 
Ap(Tm)/pN ~ 10-2 for YBaCuO and Ap(Tm)/pN ~ 
10"^ for BISCO, which should be compared with the 
jumps observed in experiments of 0.1 Q and 10~'' |19| ] 
respectively. 

An alternative estimate of the relative jump at melt- 
ing can be can be computed from the expression Px = 
exp —Ux/ksTm, and using results from experimental fits 
pO| and numerical simulations pit] that indicate Ux ^ 
kBTm- For the observed Ap(Tm)/pN — 0.1 we obtain 
Ux = '2.3kBTm, whereas the numerical simulations give 
Ux = 7-5kBTm- Even though there is a numerical dis- 
crepancy, the fact that the crossing energy is independent 
of field at the melting temperature constitutes, in our 
view, an additional element supporting the applicability 
of Verlet's criterion in the vortex state. 

Our theory asumes a continuum approximation vor the 
vortex degrees of freedom, and will be valid as long as Iz is 
larger than the interplane separation s. When Iz becomes 
of the order of s one attains the so called decoupled or 
"superentangled" regime, in which our approximations 
are no longer valid. Therefore, we expect deviations from 
the universal jump in the resistance for magnetic fields 
larger than Bx2 = g4>o/ skBT. 

An experimental proof of the universalities in S{q) 
could be obtained from neutron scattering experiments 
|18[ performed at the melting temperature, or probably 
by decorating a sample rapidly quenched from the melt- 
ing point. 

In summary, we have presented an explanation for the 
universal jump in the resistivity at the melting transition 
of vortices in high temperature superconductors. The 
theory is based on universalities at the melting transition 
as reflected in the structure factor, together with vortex 
cutting dominating the viscosity of vortices. 
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